Deformation Energy Minima at Finite Mass Asymmetry 
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A very general saddle point nuclear shape may be found as a solution of an integro-differential 
equation without giving apriori any shape parametrization. By introducing phenomenological shell 
corrections one obtains minima of deformation energy for binary fission of parent nuclei at a finite 
(non-zero) mass asymmetry. Results are presented for reflection asymmetric saddle point shapes of 
thorium and uranium even-mass isotopes with A — 226-238 and A = 230-238 respectively. 

PACS numbers: 24.75.-|-i, 25.85.Ca, 27.90.-|-b 



INTRODUCTION 

One of the earliest features observed in nuclear fission 
was the preference for breakup into two fragments of un- 
equal mass . The asymmetric distribution of fragment 
masses from the spontaneous or low excitation energy 
induced fission was a longstanding puzzle of the theory 
In the framework of the liquid drop model (LDM), 
the mass distribution of fission fragments is symmet- 
ric. By adding shell corrections within Strutinsky's Q 
macroscopic-microscopic method, it was shown that 
the outer barrier for asymmetric shapes is lower than for 
symmetric ones. In this way it was possible to explain 
qualitatively the fission asymmetry. Significant progress 
was achieved with the development of the asymmetric 
two center shell model by the Frankfurt school jaH- The 
fragmentation theory was succesful in describing both re- 
gions of low and high mass asymmetry. 

The shapes during a fission process from one parent nu- 
cleus to the final fragments, have been intensively studied 
either statically (looking for the minimimum of poten- 
tial energy) [a, Isl , or dynamically (by choosing a path 
with the smallest value of action integral) |lOl lllj . Par- 
ticularly important points on a potential ene rgy surface 
are those corresponding to the ground state [l3|, saddle 
point(s) ^] and scission point (l^l- In a static ap- 
proach, the equilibrium nuclear shapes are usually deter- 
mined by minimizing the energy functional on a certain 
class of trial functions representing the surface equation. 
The required number of independent shape parameters 
may be as high as nine values 0; they are discussed in 
ref. 0. 

The purpose of this paper is to present a method allow- 
ing to obtain a very general reflection asymmetric sad- 
dle point shape as a solution of an integro-differential 
equation without a shape parametrization apriori intro- 
duced. This equation was derived by minimizing the 
potential energy with constraints (constant volume and 
given deformation parameter) . The method 9, 15] allows 
to obtain straightforwardly the axially symmetric surface 
shape for which the liquid drop energy, Eljjm = Eg+Ec, 



is minimum. By adding the shell correction SE to the 
LDM deformation energy, E^ef = Eldm + SE, we can 
obtain minima at a finite value of the mass asymmetry 
parameter. Phenomenological shell corrections are used. 
Results for binary fission of parent nuclei 226-238rpj^ 
230-238-(j g^j.g presented. 



DEFORMATION ENERGY 

One assumes cylindrical symmetry. The deformation 
parameter a is defined as the distance between centers of 
mass of the fragments lying at the left hand side and right 
hand side of the plane z = 0, respectively: a = | z£ | -1- 1 z|j | . 
This definition allows to reach all intermediate stages of 
deformation from one parent nucleus to two fragments 
by a continuous variation of its value. The position 
of separation plane, z = 0, is given by the condition: 
{dp/dz)^^Q = 0. 

We are looking for a nuclear surface equation p — p{z) 
in cylindrical coordinates, which minimizes the potential 
energy of deformation with two constraints: volume con- 
servation, and given deformation parameter, a, assumed 
to be an adiabatic variable. 

Relative to a spherical shape, the deformation energy 
is defined by 

Ed,}{a)~E'' ^ E'l[Bs-l + 2X{Bc-l)]+5E-5E° (1) 

where E°^ = a,(l - kP)A^'^ and E^^ = acZ^A-^/^ 
are energies corresponding to spherical shape and / — 
{N — Z)/A. The relative surface and Coulomb ener- 
gies Bs = Es/E°, Be = Ec/E^ and the shell cor- 
rection SE{a) are functions of the nuclear shape. The 
dependence on the neutron and proton numbers is con- 
tained in E^, the fissility parameter X = Eq/{2E^) = 
[3Z2eV(5i?o)]/2K(l - k/2)A2/3], as well as in shell cor- 
rection energy 6E'^ of the spherical nucleus. From a fit 
to experimental data on nuclear masses, quadrupole mo- 
ments, and fission barriers, the following values of the 
parameters have been obtained ^^: — 17.9439 MeV, 
K = 1.7826, ac = 3eV(5ro) = 0.7053 MeV. The radius 
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of spherical nucleus is i?o = roA^^^ with rg = 1.2249 fm, 
and = 1.44 MeV-fm is the square of electron charge. 
The shape-dependent dimensionless surface term is pro- 
portional to the surface area, and the expression of 
Coulomb energy is a double integral. Both are com- 
puted by Gauss-Legendre numerical quadratures. 

To the LDM energy we add a phenomenological shell 
correction SE{a) by using a formula adapted after Ref- 
erence . At a given deformation we calculate the vol- 
umes of fragments and the corresponding numbers of nu- 
cleons Zi{a), Ni{a) [i = 1,2), proportional to the vol- 
ume of each fragment. Then we add for each fragment 
the contribution of protons and neutrons 

5E{a)=Y, SE, (a) = i^Epr (a) + SE^^ (a)] (2) 

i i 

given by 



SEp, = Cs{Zi)- 6Em = Cs{Ni) 



where 

s{Z) ^ Z-^/^F{Z) - cZ^/^ 
and a similar equation for s{N). 



(3) 



(4) 



F{n) 



5/3 
i-1 



(5) 

where n G (Ni-i, Ni) is the current number of protons 
(Z) or neutrons (A'') and Ni-i, Ni are the nearest magic 
numbers. The parameters c = 0.2, C = 6.2 MeV were 
determined by fit to experimental masses and deforma- 
tions. The dependence on deformation a 0| is given 
by 



SEia)^^i^Y.[siN,) 



Li{a) 
Ri 



(6) 



where Li(a) are the lenghths of fragments along the sym- 
metry axis. During the deformation process, the varia- 
tion of separation distance between centers induces the 
variation of the geometrical quantities and of the corre- 
sponding nucleon numbers. Each time a proton or neu- 
tron number reaches a magic value, the correction energy 
passes through a minimum, and it has a maximum at 
midshell. 



INTEGRO-DIFFERENTIAL EQUATION 

As mentioned above, the nuclear surface equation with 
axial symmetry around z axis is expressed as p = p{z) 
in cylindrical coordinates. We use the following rela- 
tionships for the principal radii of curvature TZi = rp. 



7^2 ~^ = —p"/t^, in which = l+p'"^. In order to mini- 
mize the deformation energy the surface equation should 
be a solution of the following equation 

pp"-p'2-[Ai + A2|z| + mXVs{z,p)]p{l + p'^f^-l = 

(7) 

where p' = dp/dz, p" = d^p/dz'^, and Vs is the Coulomb 
potential on the nuclear surface. In this equation Ai and 
A2 are Lagrange multipliers corresponding to the con- 
straints of volume conservation (or given mass asymme- 
try if the voulume is conserved in each "half" of the nu- 
cleus) and a determined value of the deformation param- 
eter a. All lengths are given in units of i?o, Coulomb 
potential in units of Ze/Ro, and energy in units of the 
surface energy E°. One can calculate for every value 
of a the deformation energy Edefio). The particular 
value Us for which dEdefip-s) I da = corresponds to the 
extremum, i.e. the shape function describes the saddle 
point (or the ground state). The associated surface equa- 
tion gives the unconditional extremum of the energy and 
corresponds to the fission barrier. The other surfaces (for 
a ^ Us) are extrema only with condition a = constant. 

The Coulomb potential on the surface depends on the 
function p{z)^ hence eq 7 is an integro-differential one, 
as Vs is expressed by an integral on the nuclear volume. 
The integration method used to solve eq 7 is based on 
the weak dependence of Coulomb energy on the nuclear 
shape. It is invariant under subtraction from Vs of a lin- 
ear function because Ai and A2 are arbitrary constants. 
The extremal surface depends on the quantity with which 
the Coulomb potential on the nuclear surface differs from 
the function Ai + A2|z|, where the constants Ai, A2 could 
be chosen in a way to minimize this difference. In the 
next iterration one uses the solution p{z) previously de- 
termined. 

The following boundary conditions have to be fulfilled 



p{zi) = p{z2) = 



lim dp{z)/dz 



00 



lim dp(z)/dz 



(8) 
(9) 



where zi and Z2 are the intercepts with z axis at the two 
tips. Equations 9 are called transversality conditions. 
For reflection symmetric shapes Z\ = —Z2 = —Zp, hence 
one can consider only positive values of z in the range 
(0,Zp). In order to get rid of singularities in eq 9 it is 
convenient to introduce a new function u{v) instead of 

u{v) = A^p'iziv)) 



where 



z{v) 



v/A 



By substituting into eq 7 one has 
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5 XV, 
2A 



4A 4A^ 



(4u + w'2)3/2] = 



(10) 



(11) 



(12) 
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Then wc introduce a linear function of v by adding and 
subtracting a + bv to 5XVs/ and defining Vsd as de- 
viation of Coulomb potential at the nuclear surface from 
a linear function of v 



Vsd = [5X/{2A)]Vs-a-vb 



(13) 



The linear term may be considered an external potential 

of deformation with a — [5X/{2A)]Vs{v = 0) and b = 
{[5X / {2A)]Va{v = Vp) - a} /vp leading to 

u"-2- l.[u'-+[(h^+a)+v{b-^) 

+Kd](4w + u'2)3/2}=0 (14) 

Here we have new constants A and Zp related to eq 10, 
besides the previous ones Ai and A2. Nevertheless the 
solution is not dependent on each parameter; important 
are the linear coefficients in v of the binomial term within 
parantheses. By equating with 1 the coefficient of v, one 
can establish the following link between parameter A and 
the Lagrange multiplier A2 



A^ = \2/^{b-l) 



(15) 



In this way u{v) should be determined by equation 

u" = 2 + i [u'^ + (v-d + V,d) (4« + u'^f'^] (16) 
containing a single parameter d. At the limit 

u(0) =0, u'(0) = l/d (17) 
and eq 9 is satisfied if Zp = Vp/A is obtained from 

u'iVpn) = (18) 

The subscript n was introduced as a consequence of the 
fact that the number of points Vpn (depending on d and 
other parameters), satisfying eq 18 is larger than unity. 
In order to solve eq 16 one starts with given values of 
parameters d and n. Different classes of shapes solutions 
of eq. 16 are obtained by taking various values of n: for 
n = 2 there is one neck (binary fission), n = 3 gives 
two necks (ternary fission), etc. For reflection symmetric 
shapes d^ — cJ_r and = rif/. Although the parameter 
A is not present in this equation we have to know it in 
order to determine the shape function from eq 10. Prom 
the volume conservation one has 



u{v)dv 



1/3 



(19) 



After solving the integro-differential equation one can 
find the deformation parameter a = 2;£ + z^, where 



= f \Ap\^)dz/ r p\z)dz 



(20) 



zI = Ia-^ £''ivp-v)uiv)dv (21) 

depends on d. From a{d), one can obtain the inverse 
function d = d{a). 

For reflection asymmetrical shapes wc need to intro- 
duce another constraint: the asymmetry parameter, 77, 
defined by 



77 : 



Ml -Mr ^ A1-A2 
Ml + Mr~ Ai+ A2 



(22) 



should remain constant during variation of the shape 
function p{z). Consequently eq 16 should be written dif- 
ferently for left hand side and right hand side. Now dL 
is different from da, and Al 7^ Ar. They have to fulfil 
matching conditions Pl{z = 0) = pr{z = 0) hence 



«L ^K)Mi = Uii"{vp)/AR 



- .1/2. 



(23) 



The second derivative p"{z) can have a discontinuity in 

z = Q \f dL ^ dR. The parameters Al and Ar are 
expressed in terms of 77, if we write eq 22 as 



2-K , 

Y 



ML = ^(l + ??)=7r^Z^ 



/ UL{v)dv 
Jo 



(24) 



2tt , 

T 



MR = =^{l-rj)=7rAR'' 



UR{v)dv (25) 



We assume that Ml+Mr is equal to the mass of a sphere 
with R — 1. From eqs 24, 25 we obtain 



^i = (l + r7)-i/M 



LO 



Ar = {1- r/)-i/MKo 
where we introduced notations similar to eq 19 



A 



LOiRO) 



2 Jo 



UL(R){v)dv 



1/3 



(26) 



(27) 



(28) 



The equation 16 is solved by successive approxima- 
tions. In every iteration one uses the 2nd order Runge- 
Kutta numerical method with constant integration step. 
The initial value u"{v = 0) can be found straightfor- 
wardly by removing the indetermination in the point 
u = 0, ?i"(0) ^ -2 + {1 - b + g)/2cF, where g = 
[5X/{2A)][dVs{v)/dv]y=o. The equation is integrated up 
to the point v = Vpn, in which the first derivative u'{vpn) 
vanishes. 



MASS ASYMMETRY 

The variations of the saddle point energy with the mass 

asymmetry parameter dL — dR (which is almost linear 
function of the mass asymmetry rj) for some even-mass 
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isotopes of Th and U are plotted in figures 1 and 2. The 
minima of the saddle-point energy occur at nonzero mass 
asymmetry parameters dj^ ~ dn between about 0.04 and 
0.085 for these nuclei. When the mass number of an 
isotope increases, the value of the mass asymmetry cor- 
responding to the minimum of the SP energy decreases. 
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FIG. 1: Saddle-point deformation energy versus mass asym- 
metry parameter for the binary fission of some even-mass Th 
isotopes in the presence of shell corrections. 
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FIG. 2: Saddle-point deformation energy versus mass asym- 
metry parameter for the binary fission of some even-mass U 
isotopes in the presence of shell corrections. 

From the saddle point energies Es p of every nucleus we 
subtract the minimum £'™p", as in the figures 1 and 2. 
The equilibrium conditions are changed when the shell 
effects are taken into account. The minimum of the Esp 
is produced by the negative values of the shell corrections 
5E - 6E°. 

In conclusion, the liquid drop model saddle point 
shapes and energy barrier heights are well reproduced by 
the present method. By adding shell corrections to the 
LDM energy we succeded to obtain the minima shown in 
Figures 1 and 2 at a finite value of mass asymmetry for 
the binary fission of 226-238rp]^ ^^^^ 230-238^ nuclei. 

As mentioned by Wilkins et al. calculations of 

PES for fissioning nuclei "qualitatively account for an 
asymmetric division of mass" . From the qualitative point 



of view the results displayed in Figures 1 and 2 proove 
the capability of the method to deal with fission mass 
and charge asymmetry. 

The experimentally determined mass number of the 
most probable heavy fragment "^o] for U isotopes ranges 
from 134 to 140. The corresponding values at the dis- 
played minima in Figures 1 and 2 are very close to 125, 
which means a discrepancy between 6.7 % and 10.7 % 
for Ah ■ The inaccuracy in reproducing the experimental 
mass asymmetry is due to the contribution of the phe- 
nomenological shell corrections. In the absence of shell 
corrections the pure liquid drop model (LDM) refiection- 
symmetric saddle point shapes [8] are reproduced, and 
the barrier height increseas with an increased mass asym- 
metry. When the shell corrections are taken into account 
the LDM part behaves in the same manner (larger values 
at non-zero mass asymmetry). Only the contribution of 
shell effects can produce a minimum of the barrier height 
at a finite value of the mass asymmetry. One may hope 
to obtain a better agreement with experimental data by 
using a more realistic shell correction model. 
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